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Three-dimensional Weyl fermions are found to emerge from simple cubic lattices with staggered 
fluxes. The mechanism is to gap the quadratic band touching by time-reversal-symmetry-breaking 
hoppings. The system exhibits rich phase diagrams where the number of Weyl fermions and their 
topological charge are tunable via the plaquette fluxes. The Weyl semimetal state is shown to be 
the intermediate phase between non-topological semimetal and quantum anomalous Hall insulator. 
The transitions between those phases can be understood through the evolution of the Weyl points 
as Berry flux insertion processes. As the Weyl points move and split (or merge) through tuning 
the plaquette fluxes, the Fermi arcs and surface states undergo significant manipulation. We also 
propose a possible scheme to realize the model in ultracold fermions in optical lattices with artificial 
gauge fields. 
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I. INTRODUCTION 

Massless Dirac fermions possess chiral symmetry and 
can be classified by their chirality 7. For example, in 
three-dimension (3D), if the Hamiltonian is H = vf<j ■ q 
where cr is the Pauli matrices vector and q is the wavevec- 
tor, then the chirality is 7 = sgn(«i?). Such fermions 
(also called chiral or Weyl fermions) possess several pe- 
culiar behaviors such as Adler-Bell-Jackiw anomaly p}. 
Chiral fermions in two-dimension (2D) have been found 
in real condensed matter systems since the discovery 
of grapheneQ. Many novel electronic properties of 
grapheneQ: Klein tunneling, the peculiar integer quan- 
tum Hall effect, the transport properties such as the 
conductivity minimum, the weak (anti) localization, and 
edge states, originate from the chiral fermion nature. 
Recent studies also found the realizations of 2D chi- 
ral fermions in ultracold atomic gases optical lattices [1]. 
A generic route to chiral fermions is to search systems 
with two-band touching^. Around the band touch- 
ing node, the Hamiltonian can generally be written as 
H = h(k) • cr where Pauli matrices cr act on the space 
spanned by the two bands and |h| — > at the node. 
For example, in graphene, around one of the node K, 
h(k) = vp(k — K). Thus chiral fermions emerge as 
topological defects (vortices) in k-space. They are gen- 
erally classified by their vortex winding number N w . In 
the special cases when N w = ±1, the winding number 
gives the chirality 7 = N w . In general cases, N w can 
be any integer. The total winding number in the sys- 
tem must be conserved under adiabatic transformation. 
The winding number can be defined by the Berry phase 
carried by the node, N w = ~ / r dk ■ (*(k)|iV k |*(k)), 
where T is a contour enclosing the node and ^(k) is the 
single- valued and continuous wavefunction of the eigen- 
states with i?|*(k)) = ±|h(k)||*(k)). By breaking the 



time-reversal symmetry, such two-band touching can be 
gapped, leading to a quantum anomalous Hall (QAH) 
insulator with Chern number C = ±4iV w p|. 

3D Weyl fermions are more robust: they can not even 
be gapped by time-reversal symmetry breaking. In fact 
they can only be annihilated in pairs with the total wind- 
ing number conserved under adiabatic transformations. 
The winding number of 3D Weyl fermion is defined asQ 

N w = ^e vbp § s dS p h- (d ku hxd ks n\ (u, S,p = x,y,z and 

s being the Levi-Civita tensor) where S is a surface en- 
closing the band touching node (Weyl point), dS p is the 
surface area elements along p direction, and n = h/|h|. 
3D Weyl fermions arc monopoles of the Berry-phase 
gauge fields where the monopole charge is the topological 
charge A^Q. As a consequence, there is a step change 
in the Hall conductivity, e.g., o~ xy (k z ) [and the Chern 
number C xy (k z ) 



■%z<7 xy (k z )} as function of k z 



<Jx V {k z ) = sgn(fc 2 - k z ) 
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due to the monopole at k z , while ... denotes other con- 
tributions to the Hall conductivity (Chern number). Ac- 
cording to the Nielsen- Ninomiya theorem^, Weyl points 
must appear in pairs with opposite N w in a lattice sys- 
tem. 

Weyl semimetal (system with 3D Weyl fermions) pos- 
sesses very special properties such as chiral surface states 
with open Fermi surfaces (Fermi arcs) terminating at 
the projection of the Weyl points and thickness depen- 
dent quantized anomalous Hall conductivity in thin films, 
which are first found in the studies of 3 He-A0, H|. Af- 
ter that Murakami showed that a Weyl semimetal phase 
can appear as intermediate phase between normal insula- 
tor and topological insulator (T^). Recent studies demon- 
strate that Weyl semimetal can also be realized in other 
condense matter systems: llj some pyrochlore iridates 
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(such as Y 2 Ir 2 07)[l2|, super lattices made of topolog- 
ical insulator and non-to polo gical insulator thin films 
with broken time- reversal jl3t Il4j or inversion fl5j sym- 
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metry, the ferromagnetic compound HgCr9Seii[l6l [13], 
and bulk magnetically doped Bi2Se 3 [l8j]. It is also found 
that in some situations the number and type of the 
Weyl points are determined and protected by the lat- 
tice symmetries (T7I [l9| . A topological nodal semimetals 
where the bulk spectrum exists nodal lines are also pro- 
posed and studied in Ref. [3| ■ There are also some lattice 
models where Weyl fermions are found poj. 

In this work we show that Weyl fermions can emerge 
from simple cubic lattices with staggered fluxes through 
plaquettes [see Fig. 1]. Differing from previous stud- 
ies where the Weyl semimetal phases emerge due to in- 
verted bands with spin-orbit interactions flf l 15 ^ 17 1 or 
gaping Dirac cones by various mass terms[l3l Il4l lsj] . 
here the mechanism is to gap the quadratic band touch- 
ing via time-reversal-symmetry-breaking hoppings due to 
the staggered fluxes. A direct distinction is that such sce- 
nario does not need to invoke spin degeneracy breaking 
(i.e., spin-orbit coupling). The system exhibits rich phase 
diagrams where the number of Weyl fermions and their 
topological charge are tunable via the fluxes per plaque- 
ttes. The Weyl semimetal state is demonstrated as the 
intermediate phase between non-topological semimetal 
and quantum anomalous Hall insulator. The transitions 
between those phases can be understood via the evo- 
lution of the Weyl points as Berry flux insertion pro- 
cesses [see Sec. IV and Summary section]. As the Weyl 
points move and split (or merge) through tuning the pla- 
quette fluxes, the Fermi arcs and surface states undergo 
significant change, because the Fermi arcs act as Dirac 
strings which have to connect the monopoles with oppo- 
site charges. Finally, we propose a possible scheme to 
realize such model in ultracold fcrmionic gases in optical 
lattices with artificial gauge fields. 



II. LATTICE AND HAMILTONIAN 

We consider a simple cubic lattice system, which can 
be viewed as stacking of layers of 2D lattices with 
checkerboard-patterned staggered fluxes [see Fig. 1] [2TI - 
[24[. The 2D checkerboard lattice is designed in such a 
way that the hopping between the nearest-neighbored A- 
type sites along the x and y directions are t x and t y 
respectively, whereas those for B-type sites are t y and 
t x respectively [see Fig. 1(a)]. The hopping from A-type 
site to the nearest B-type ones are t^t ^ 1 as indicated 
in Fig. 1(a). The flux per plaquette is ±$1 = ±40i 
[Fig. 1(a)]. Recently Sun et al. proposed a realistic 
optical lattice system to realize such model 21]. Be- 
sides it can be realized in ultracold fermions in op- 
tical lattices and in condensed matter systems with 
artificial [25l-[29j or "emergent" [3(| gauge fields. The 
Hamiltonian for each layer is |2lTl24l | H2B = ho(k)ao + 
hi(k) • cr where <x is the Pauli matrix vector acting 
on A/B (pseudo-spin up/down) site space and (7 is 
the 2x2 identity matrix. As it is no need to in- 
voke true-spin splitting, we keep the true-spin states 



as degenerate. h Q (k) = 2t Q (cosk x + cosk y ), h lz {k) = 
2t\ (cos k x — cos k y ) , hi x (k) = 4t 2 cos 4>i cos cos , and 

/iiy(k) = 4t2 sin 4>i sin 4^ sin ^ where to = (t x + t y )/2, 
t\ = (t x — t y )/2 [Note that throughout this paper, we 
set the lattice constant a = 1. All the wavevectors are 
then in the units of 1/a = 1.]. The two bands touches 
quadratically at (jr, tt) when sin cii = 0. This band touch- 
ing is nontrivial as it carries a nonzero winding number 
N w = 2sgn(iii2 cos <fii) = ±2. At finite sincii (broken 
time-reversal symmetry) the quadratic band touching is 
gapped and the system becomes a QAH insulator with 
Chern number C = hsgjitit sin <b\ )N W [221 [23j . 
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FIG. 1. (Color online) Lattice structure in (a) (001) plane 
in an odd layer and (b)(110) plane. Red-dots and blue- 
squares represent type A and B sites respectively. Hoppings 
along (opposite) the green-solid and blue-dotted arrows are 
t 2 e" i01 and t ± e"^ 2 (t 2 e <01 and t^ 2 ) separately. Those 
along red-dashed and blue-dotted lines are t v and t x respec- 
tively. ±$1 and ±$2 are the plaquette fluxes. Phase diagram 
at (c) \t 2 /t±\ = 0.95 and (d) \t 2 /t±\ = 1.3. Yellow region: 
Weyl semimetal; green region: QAH state with Chern num- 
ber C xy (k z ) = — 2sgn(tit2t±) for all k z ; brown region: QAH 
state with C xy (k z ) — 2sgn(tit2t±) for all k z . 

The 3D lattice is a stacking of the 2D layers in such 
a way that different types of sites are on top of each 
other [see Fig. 1(b)]. The hopping between those sites is 
t±e ± ^ 2 . The staggered flux per plaquette in the (110) 
plane is ±$ 2 = ±2((/>i + (/>2) [in (110) plane, the plaquette 
flux is ±2((/>i — (j>2)\- The Hamiltonian of the system is 
then 



H = Mk)<7 + h(k) • o-]r 

+ g{k z ){cos(t)2crx + sin 4>2cr y )n 



(2) 



Here a and t matrices act on the A/B sites and odd/even 
layers respectively. g(k z ) = 2tj_ cos 4f . The Hamiltonian 
can be block diagonalized directly, giving 



H± = h {k)a + h±(k) ■ cr 



(3) 



3 



in the two blocks where 



k x ky 



h± x (k) = 4i 2 cos0i cos — cos ± g{k z ) cos 02, 
k k 

h±y(k) = 4^2 sin 0i sin sin ± g(k z ) sin 02, 
h± z (k) = 2ti(cosk x — cos k y ). 



(4) 



For the sake of easier formulation, we set the first Bril- 
louin zone as k x <E [0, 2tt], k y G [0, 2tt] and k z € [0, 2tt]. 
Finally, the system possesses a series inversion symme- 
tries. It is invariant under the following inversion trans- 
formations: (i) z — > —z, (ii) (x,y) — > (— x, — y), and (hi) 
(x,y, z) — > (— x, —y, —z). At ho = 0, the system also has 
particle-hole symmetry 



III. PHASE DIAGRAMS 
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The energy spectrum of the system is 
£7 a ±(k) = h (k)±|h a (k) 



(5) 



with a = ±. The equations for the nodes (Weyl points) 
are h ax = h ay = h az = 0. One finds that away from 
0i, or 2 = 0, ±7r/2, 7r [see Fig. 1 and discussions below], 
there are four Weyl points K + ^j = (f3k^.,f3k y , k z ) (/3 = ±) 
in the + block and K-^ = (/3fc£, (3k y , 2tt — k z ) in the — 
block where 



k y = 7r — 2 arctan - 



tan 0i 



tan 02 



, r 7 r • t Z' tan 01 \ 

k x = V k y with J7 = Bg a^_j > 

tan 0i 



k z = 2 arccos 



— 277^2 sin 0i 



tan 02 



t_L sin 02 V 
Note that the above equation holds only when 



2t 2 sin 0i / 
I 1 

t_L sin 02 V 



tan 0i 



tan 02 



< 1, 



(0) 



(7) 



which sets the phase boundaries of the Weyl semimctal. 
It is required that \ho(k)~- ho(K a p)\ < |h Q (k)| so that the 
system is an insulator wherever away from the nodes K Q) g 
(the definition of semimetal). Given this, the ho(k)ao 
term is irrelevant for the physics to be discussed and we 
hence take ho = hereafter. 

In Fig. 1(c) and 1(d) we plot the phase diagram at 
two different ^A-lI- There are three phases: the Weyl 
semimetal, the QAH, and the non-topological semimetal. 
The last one refers to semimctals which have only nodes 
with N w = 0. Such nodes are accidental band degeneracy 
points, which can be gapped by infinitesimal band mix- 
ing without breaking any symmetries. The QAH phase 
actually consists of two topologically distinct phases with 
opposite Chern number [C xy (k z ) — ±2sgn(tii 2 i jj f° r 
all k z ]. Those two phases are separated by the Weyl 



FIG. 2. (Color online) (a) and (c) Illustration of the Weyl 
points when there are four in k y -k z plane. The arrows in- 
dicates the 2-direction movement of the Weyl points as 0i 
increases. Blue-squares (red-dots) denote Weyl points with 
N w = —1 (N w = 1). The parameters are ti = 1, i 2 = 2, 
t± = 2, and 02 = — 0.25-7T. cj>\ — — O.ln in (a), and cj>i = Q.ln 
in (c). (b) and (d) Chern number C xy (k z ) as function of k z 
for the cases in (a) and (c) respectively. 



semimetal and non-topological semimetal phases. In 
fact, the Weyl semimetals can be viewed as the inter- 
mediate phases between the non-topological semimetal 
(or insulator) phase and the QAH phase. In the for- 
mer the Chern number is always zero, whereas in the 
latter it is always nonzero. In Weyl semimctals, there 
are two regions in the Brillouin zone: the Chern number 
is zero in one region while nonzero in another. Specifi- 
cally in the current model C xy (k z ) is nonzero only when 
k z £ (kf in ,kf ax ) with kf iD = Mm(k c z ,2n - kf) and 
kf ax = Max(fc£, 2n-k c z ). The QAH and non-topological 
semimetal phases can be viewed as the limit fc™ ln — > 
and fc" lln — > 7T respectively. There are actually two Weyl 
semimetal phases in which the Chern number C xy (k z ) 
is opposite at k z £ (fc" lm , fc™ ax ). Those two states are 
the intermediate states between the two QAH states 
(with opposite Chern number) and the non-topological 
semimetal state separately. The latter lies at the lines 
0i = 0,±7r/2, and n in the phase diagrams, where 
k" un = /c" lax = 7T. We will show in next section that 
the evolution of the ground state and quantum phase 
transitions between those phases can be understood via 
the evolution of the Weyl points as Berry flux insertion 
processes. 

We find that as I^A-i-l increases, the area of the Weyl 
semimetal phase in the phase diagram as functions of 
0i and 02 shrinks. At I^A-iJ — * °°> such area be- 
comes zero. When ^A-i-l — 1/2 the system is al- 
ways in the Weyl semimetal phase except at the spe- 
cial lines 0i = 0, ±7r/2,7r. Hence the region of the Weyl 
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scmimetal phases can be tuned by the ratio |*2/*x|- It is 
noted that the phase diagram exhibits some angle-shaped 
structures around the special points (fa, fa) = (fa, ±fa) 
with fa = 0,±7r/2,7r as well as (0, tt) and (n, 0) where 
Eq. ([6]) becomes ill-defined. The structure around, say, 
(0,0), can be understood via the following analysis. For 
(fa, fa) close to (0,0), along the line fa = £fa, k z = 
2 arccos[— 2*2/*x(|£| + 1)]. Hence the Weyl semimetal 
phase is at |£| > 1 2*2 / 1 j_ | — 1, which is angle-shaped. 
When |*2 Ax I < 1/2) the system is in the Weyl semimetal 
phase for all parameters (fa, fa) around (0,0). 

It is benefit to point out some special regions in the 
phase diagram. First, when fa = ±7r/2, there are only 
two Weyl points: K+ = (tt, tt, kf) in the + block and 
K_ = (tt, tt, 2tt — k z ) in the — block. Those Weyl points 
are quadratic band touchings with winding number N w = 
±2. Similarly, when fa — 0, tt, there are two quadratic- 



band-touching Weyl points: K + = (0,0, k z ) in the + 
block and K_ = (0, 0,27r— k c z ) in the — block. Besides, 
as pointed out before at the lines fa = 0,±7r/2,7r the 
system is a non-topological semimetal. 



IV. EVOLUTION OF WEYL POINTS 

Away from the above regions, there are four Weyl 
points, around which fermions are described by the Weyl 
Hamiltonian 

H aP (k) = (T-v a p-q+0(q 2 ) (8) 

where q = k K Q ^. The velocity tensors of the Dirac 
cones are 



■/3r]t2 cos fa sin ky — /3*2 cos fa sin k c y —at 



/3*2 sin 



>i sin k y 
2/3r/iisin k c y 



(3rjt2 sin fa sin ky 
2/3ti sin k y 



. cos fa sin 
—at± sin fa sin -f- 




(9) 



The winding number is the sign of the determinator of 
the velocity tensor N w (a(3) = sgn[dct(u Q , /3 )]. One finds 

N w (af3) = -asgn[*i*2*_L s'm(fa + rjfa)}. (10) 

Note that Weyl points in the same block a have the same 
topological charge (independent of /3), whereas Weyl 
points in different block have opposite topological charge. 
The positions and motions of the four Weyl points are il- 
lustrated in Fig. 2 for a specific case where N w = —a, 
ky = —k c x < tt, and k c z > tt. According to Eq. ([T]), the 
Chcrn number C xy (k z ) varies with k z . The Chern num- 
ber changes only when the gap is closed and re-opened, 
i.e., when k z passes through k c z and 2n — k z . The Chcrn 
number is calculated as 

C xy = ^ — / dk x / dk y n a ■ (d k:i n a x d ky n a ) , 

a=± 47F •'O ^0 

(11) 

where n a = h Q /|h a |. We find that C xy (k z ) = 
2sgn[*i* 2 * ± sin(2(/)i)], when k z G (fcf in ,fcf ax ), otherwise 
C xy (k z ) = 0. The relation between the quantum phase 
transitions and the evolution of the Weyl points is de- 
picted in Fig. 2 where we consider a situation with 
fa = — 0.25-7T and fa varying from — 0.17T to 0.17T. At 
4>i = — 0.1-7T, the Weyl points with N w < are in the 
k z < tt region whereas those with N w > are in the 
k z > tt region [Fig. 2(a)]. Consequently, the Chern num- 
ber C xy (k z ) is negative at k z £ (fcf in ,fc™ ax ). As fa in- 
creases the Weyl points with N w < and those with 
N w > moves in opposite directions along k z and be- 
comes closer. At fa = 0, Weyl points with opposite 
N w coincide and the system becomes a non-topological 



semimetal. After that their positions shift as fa in- 
creases [Fig. 2 (c)]: the Weyl points with N w < in 
the k z > tt region whereas those with N w > move 
to the k z < tt. As a consequence the Chcrn number 
C xy (k z ) changes sign at k z € (fc™ ln , fc™ ax ). Further vari- 
ation of fa will enlarge the region and finally the system 
becomes a QAH insulator after fc™ ln — > where pairs of 
monopoles with opposite topological charge (N w ) merge 
and annihilate each other. During those processes quan- 
tized Berry fluxes are inserted into each k x -k y plane with 
fixed k z , whenever the monopoles move across it, as there 
are quantized Berry fluxes flow between monopoles with 
opposite charge. 

In addition, as fa — ¥ 0,±7r/2, and tt, the two Dirac 
cones with same k z merge together and form a quadratic 
band touching as they have the same winding number. 
The k z dependence of C xy (k z ) in this situation is similar 
to the previous one. However, as we will show later, the 
Fermi arcs and surface states in those two situations are 
significantly different. 



V. EVOLUTION OF FERMI ARCS AND 
SURFACE STATES 

The move and merge (or split) of the Weyl points have 
profound effects on the topologically protected surface 
states. In particular, we demonstrate the Fermi arcs on 
the surface for two situations where there are (i) four 
Dirac cones and (ii) two quadratic band touchings in 
Fig. 3 (a) and (b) respectively. The color represents the 
spectral function A(E) = ^-ln\G r (E) [G r (E) is the re- 
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FIG. 3. (Color online) Spectral function A(E) of the surface 
states at zero energy (E = 0) in the case with (a) four Dirac 
cones and (b) two quadratic band touchings. Spectra of sur- 
face and bulk states as function of (c) k z (at k y = 0.8n) and 
(d) k y (at k z = 7r) respectively with the same parameters as 
in (a). The size of the system is 151 unit cell (due to the fi- 
nite size effect the bulk bands are slight gapped at the node). 
Gray region represents bulk spectra, while red (green) curves 
denote the surface spectra at the left (right) boundary [Note 
that in (c) the surface spectra at the two boundaries coincide]. 
ti = 1, ta = 2, and t± = 2. In (a), (c) and (d) (j>i — 0.1-71", 
<j>2 = — 0.257T, whereas in (b) <j>i = 0.157T, (f>2 = — 0.5-7T. Cor- 
respondingly, k y = 0.8-7T and k% = lAn in (a), (c) and (d), 
whereas k y = tt and k c z = 1.57T in (b). The results are cal- 
culated via the iterative Green function method in Ref. [3l| . 
An artificial spectral broadening 0.04 is taken for the sake of 
visibility. 



tarded Green function of the system] of the zero-energy 
(E = 0) surface states (i.e., the Fermi arcs). The surface 
is perpendicular to the x direction. It is noted that al- 
though the k z dependence of C xy (k z ) are similar in the 
two cases, the Fermi arcs are quite different. This is es- 
sentially due to the different positions of the monopolcs 
in k-space. From Fig. 3(c), one can see that the zero- 
energy surface states (Fermi arcs) merge into the bulk 
bands through the Weyl points. This figure clearly visu- 
alizes the fact that Fermi arcs are Dirac strings which 
link the monopoles and antimonopolcs of opposite "mag- 
netic" charge (N w ) in the Berry-phase gauge fields p} . 
As Fermi arcs terminate at the Weyl points, the posi- 
tions of the monopoles are essential to the shape of the 
Fermi arcs and the topologically protected surface states 
spectra. From the figure it is seen that the monopole- 
antimonopole pair with the same k y (but opposite k z ) 
are connected via the Dirac strings. This string configu- 
ration is due to the facts that (i) the system possesses the 
inversion symmetry: k z — > — k z \ (ii) the string connect- 



ing the monopolcs and anti-monopoles must pass through 
the k z = 7T plane because it exists only in the region 
k z e (kf in , kf ax ) dictated by the topology. Hence the 
string connection between monopole-antimonopole pair 
with the same k y (and opposite k z ) is protected by the 
symmetry and dictated by the topology. 
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FIG. 4. (Color online) Spectral functions of the surface states 
for two cases: (a) four Dirac cones (b) two quadratic band 
touchings. The energies of the surface states from left to 
right (in both the k y < ir and k y > ir regions) are -0.8, -0.4, 
0, 0.4, and 0.8 respectively. E.g., both the leftmost arc in 
the k y < 7r region and the leftmost arc in the k y > tt region 
are the arcs with energy -0.8. The parameters for (a) and 
(b) are the same as those in Fig. 3 (a) and 3 (b) respectively. 
An artificial spectral broadening 0.04 is taken for the sake of 
visibility. 

It is noted from Fig. 3(b) that in case (ii) the two Fermi 
arcs are actually connected together as ky = — ky, form- 
ing a closed Fermi surface. In the four Dirac cone case in 
Fig. 3(a), when ky ^ — ky, the Fermi surface is not closed. 
It seems that the Fermi arcs are also invariant under the 
inversion operation k y — > — k y . However, this is only a 
special property of zero energy surface states due to the 
particle-hole symmetry. For surface states with nonzero 
energy, the spectral function are not fc y -inversion sym- 
metric as the boundary breaks the (x, y) — > (— x, —y) in- 
version symmetry. To illustrate the surface states explic- 
itly, we plot the spectral functions of the surface states 
in Fig. 4(a) and 4(b) for the two cases that correspond 
to Fig. 3(a) and 3(b) respectively. The selected ener- 
gies of the surface states from left to right (in both the 
k y < tt and k y > tt regions) arc -0.8, -0.4, 0, 0.4, and 
0.8 respectively. For each energy there are two arcs: one 
in the k y < tt region and another in the k y > tt region. 
E.g., both the leftmost arc in the k y < tt region and the 
leftmost arc in the k y > tt region are the arcs with en- 
ergy E = — 0.8. It is clearly seen that the spectra are 
not invariant under fc^-inversion but invariant under the 
particle-hole transformation. Besides, only the arcs with 
zero energy links between the monopoles whereas other 
arcs merge into the bulk bands without going through 
the Weyl points. Finally, it is seen from Fig. 3(d) that 
the two Fermi arcs are assigned to two different chiral 
edge states at the left boundary which both have posi- 
tive group velocity along the z-direction. It is found that 
the shape of the Fermi arcs can be tuned via <f>\ and <p2 
to be inner-curved (as in Fig. 3b), out-curved, or flat [as 
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in Fig. 3(a)]. Flat bands can be interesting in the context 
of elevated transition temperature in spontaneous sym- 
metry broken 0] , as the density of states are increased. 



VI. A POSSIBLE SCHEME FOR 
EXPERIMENTAL REALIZATION 



In this section we propose a possible scheme to realize 
the model in optical lattices. The required artificial gauge 
fields are generated by the spatial variation of the laser- 
atom interaction as suggested in Ref. H3 ■ Suppose that 
there is a excited states with energy much higher than the 
energy scale where the above model are defined. Impose 
standing waves of light to induce the following coupling 
between the ground and excited states in rotating wave 
approximation, 

Hr{t) = M[cos2-k z cos it [x+y)F x +cos it {x—y)F y +C > F z \ 1 

(12) 

where F v (y — x,y,z) are the Pauli matrices acting on 
the (dressed-) ground and excited states. M and ( are 
parameters of the laser-atom coupling. It can be real- 
ized in a system with three standing wave Raman lasers 
with a detuning of MCH^]. The three laser wavevec- 
tors are (it, tt,±2tt) and (it, — ir,0). When the kinetic 
energy is small compared to the energy splitting of the 
local dressed states [the eigenstates of the local Hamil- 
tonian _ff^(r)]. The emergent gauge fields can be ob- 
tained via the Berry-phase |27j . A = ( 1 i'G(r)\i'V r \ 1 i'G( r )) 
with \^a( r )) being the ground state of the local Hamil- 



tonian Hji(r) 
Q = 1 



One can show that A 
ICI 



ieV r II where 
and II = 



cos 2 (2"7rz) cos 2 7r(o;-t-y)+cos 2 7r(x — y)+C 

Arg[cos(27rz) cos7r(a; + y) + zcos7r(x — y)]. The effective 
"magnetic field" (i.e., the Berry curvature, or the "mag- 
netic flux density") are written as B v = e" Sp dsA p , with 
v, 5, p ~ x, y, z and s being the Levi-Civita tensor. The 
plaquette flux for each plaquette is the "magnetic flux" 
through it. For example, for a plaquette in the x-y plane, 
the plaquette flux is $ — J dxdyB z , where the integral is 
limited within the plaquette. We plot the effective "mag- 
netic field" B z through the x-y plane (in odd layer) and 
£?[iTo] through the (110) plane in Fig. 5(a) and 5(b) re- 
spectively It is seen that the magnetic flux density has 
exactly the same checkerboard pattern as the plaquette 
fluxes in Fig. 1(a) and 1(b). This shows that the gauge 
fields generated by the Hamiltonian, Eq. (fT2"j). is exactly 
what is needed for the realization of the model. The 
magnetic flux density in (110) plane is zero everywhere, 
which is because the current scheme realizes the model 
with fa = fa. 

The phase along a hopping path I is given by fa = 
^dr ■ A. Combined with Figs. 1(a) and 1(b), one can 
show that the hopping phases fa and fa in Hamiltonian 
Eq. ^ are given by 



FIG. 5. (Color online) The "magnetic flux density" through 
(a) (001) plane at an odd layer (b) (110) plane, r = r ■ 
(1, 1, 0)/2 = (x + y)/2. The bright region has positive value, 
whereas the dark region has negative value. 



whereas the phases of other hoppings are zero. The func- 
tion g(C,) is plotted in Fig. 6. fa = fa € [— 7r/4, 0]. 
Due to the factor cos(27T2:) the gauge phases are inverted 
from odd layer at z = n to even layer at z = n + i 
with n being integer. This exactly realizes the lattice 
structure in Fig. 1. The amplitudes and signs of the 
rest hoppings can, in principle, be tuned via various op- 
tical lattice techniques f32|. By manipulating the ratio 
ia/^X) the system can experience various phases in the 
phase diagram: (i) At finite <fii when |t2/t±| < 1 it is a 
Weyl semimetal, (ii) otherwise it is a QAH insulator with 
C xy {k z ) = 2sgn[tit 2 t± sin(2(/>i)] for all k z \ (iii) At fa = 
the system is always a non-topological semimetal. 




fa = 02 = g(0 



(13) 



FIG. 6. The function g((). 



VII. SUMMARY 

In summary, a model of simple cubic lattice with stag- 
gered fluxes which exhibits Weyl semimetal phase is pro- 
posed and studied. The model is simple and can act as 
a prototype to study the properties of Weyl semimctals. 
Due to its simplicity, the model is potentially achiev- 
able both in ultracold fcrmions in optical lattices and 
in condensed matter systems. In particular, we propose 
a possible scheme to realize the model in optical lattice 
system. Differing from previous works, here the mech- 
anism to achieve the topological Weyl semimetal state 
is to gap the quadratic band touching by time-reversal- 
symmetry-breaking hoppings. The system exhibits rich 
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phase diagrams, where the number of Weyl fermions and 
their topological charges and positions arc tunable via the 
plaquette fluxes (hopping phases). The Weyl semimctal 
state is demonstrated to be the intermediate phase be- 
tween non-topological semimetal and quantum anoma- 
lous Hall insulator. The transitions between those phases 
can be understood via the evolution of the Weyl points 
[sec below] . As the Weyl points move and split (or merge) 
via the manipulation of the hopping phases, the Fermi 
arcs and surface states undergo significant change, as the 
Fermi arcs have to be terminated at the Weyl points. 

The relations between the non-topological insula- 
tors/semimetals, topological Weyl semimetals and QAH 
insulators demonstrated in this paper can be summa- 
rized in the following processes: (i) A non-topological 
insulator becomes non-topological semimctal via forming 
an accidental band-touching node with winding number 
N w = 0. (ii) By splitting the node into pairs of Weyl 
fermions with opposite topological charges and moving 
the positively and negatively charged Weyl fermions in 
opposite directions in the k-space, a Weyl semimctal 
state is created, (iii) When the positively and negatively 
charged Weyl points are moved by half of the recipro- 
cal lattice vector, they merge and annihilate each other 
in pairs, making the system transit into a QAH insula- 



tor. The sign of the Chern number of the QAH state de- 
pends on the direction along which the Weyl fermions are 
moved. Quantized Berry fluxes are inserted into the bulk 
states during the evolution because Weyl fermions are the 
monopoles of Berry-phase gauge fields. This picture can 
be generalized to understand the relations between the 
non-topological insulators/semimctals, the time- reversal 
symmetric topological Weyl semimetals and the 3D quan- 
tum spin Hall insulator (i.e., the time-reversal invariant 
Z2 topological insulators), where quantized non-Abelian 
Berry fluxes are inserted into the bulk states during the 
moving of the monopoles. 
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